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The authors have introduced recently a "micro canonical functional integral" which 
yields directly the density of states as a function of energy. The phase of the functional 
integral is Jacobi's action, the extrema of which are classical solutions at a given 
energy. This approach is general but is especially well suited to gravitating systems 
because for them the total energy can be fixed simply as a boundary condition on the 
gravitational field. In this paper, however, we ignore gravity and illustrate the use 
of Jacobi's action by computing the density of states for a nonrelativistic harmonic 
oscillator. 
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2 INTRODUCTION 

Jacobi's form of the action principle involves variations at fixed energy, rather than 
the variations at fixed time used in Hamilton's principle. The fixed time interval 
in Hamilton's action becomes fixed inverse temperature in the "periodic imaginary 
time" formulation, thus transforming Hamilton's action into the appropriate (imagi- 
nary) phase for a periodic path in computing the canonical partition function from a 
Feynman functional integral (Feynman and Hibbs 1965). In contrast, fixed total en- 
ergy is suitable for the micro canonical ensemble and, correspondingly, Jacobi's action 
is the phase in an expression for the density of states as a real-time "microcanonical 
functional integral" (MCFI) (Brown and York 1993b). 

We wish to characterize briefly the canonical and microcanonical pictures. (We shall 
speak only of energy and (inverse) temperature here, ignoring the other possible 
conjugate pairs of variables in order to simplify the discussion.) In the canonical 



picture, with a fixed temperature shared by all constituents of a system, there are no 
constraints on the energy. This feature simplifies combinatorial (counting) problems 
for canonical systems and leads to the factorization of the partition function for 
weakly coupled constituents. For gravitating systems in equilibrium, the temperature 
is not spatially uniform because of gravitational red and blue-shift effects. In such 
cases the relevant temperature is that determined at the boundary of the system 
(York 1986). It can therefore be specified as a boundary condition on the metric 
(York 1986, Whiting and York 1988, Braden et al. 1990) and used in conjunction with 
Hamilton's principle, which is the form of the action for gravity in which the metric 
is fixed on the boundary (Brown and York 1992, 1993a). (The metric determines 
the lapse of proper time along the boundary.) On the other hand, equilibrium in the 
canonical picture is not always stable when gravity is present, as is well known. For 
some pertinent examples, see York (1986), Whiting and York (1988), and Braden et 
al. (1990). 

With its constraint on the energy, the micro canonical picture leads to more robust 
stability properties. However, the energy constraint can complicate calculations of 
relevant statistical properties because the constituents of the system share from a 
common fixed pool of energy. For field theories, with a continuous infinity of degrees 
of freedom, the energy constraint restricts the entire phase space of the system un- 
less gravity is taken into account. For gravitating systems, as a consequence of the 
equivalence principle, the total energy including that of matter fields is an integral 
of certain derivatives of the metric over a two-surface bounding the system. There- 
fore, if we specify as a boundary condition the energy per unit two-surface area, we 
have constrained the total energy simply by a boundary condition (Brown and York 
1993a, 1993b). Thus, the canonical and microcanonical cases differ only in which of 
the conjugate variables (Brown et al. 1990), inverse temperature or energy, is speci- 
fied on the boundary. The corresponding functional integrals, for partition function 
or density of states, differ in which action gives the correct phase, Hamilton's or 
Jacobi's. 

We have recently applied this reasoning to the case of a stationary black hole (Brown 
et al. 1991a, 1991b, Brown and York 1993b). The MCFI, in a steepest descents 
approximation, shows that the logarithm of the density of states is one-quarter of 
the area of the event horizon (that is, the Bekenstein-Hawking entropy) (Brown and 
York 1993b). In the present paper we shall disregard gravity and obtain the density 
of states for a nonrelativistic harmonic oscillator. This is a relatively simple situation 
in which to recall the properties of Jacobi's action and to see the MCFI at work. 
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3 JACOBI'S ACTION 

Consider, for simplicity, a particle of mass m with a one-dimensional configuration 
space. The Lagrangian form of Jacobi's action is (Lanczos 1970, Brown and York 
1989) 

S E [x] = J dx^2m[E-V(x)] , (1) 

where V(x) is the potential energy and the energy E is a fixed constant. S E [x] is 
extremized by varying the path freely except that the end points are fixed. Now 
introduce a parameter a increasing monotonically from a' at one end of the path to 
a" at the other. Denoting dx/da by x, we can write the action as 

S E [x] = J dax^2m[E -V(x)] , (2) 

where x' = x(a') and x" = x(a") are fixed. Jacobi's action is invariant under changes 
5x induced by changes of parameterization that preserve the end-point values of a. 

For constructing the MCFI, we employ the canonical form of Jacobi's action. Because 
of the reparameterization invariance of S E [x] , the corresponding canonical Hamilto- 
nian x(dL/dx) — L vanishes identically. Furthermore, the canonical momentum 

P =^=[2m(E-V)} 1/2 (3) 

is independent of x in one dimension and, in general, does not allow one to solve 
for all the i's as functions of the p's. Indeed, from (3) we obtain the "Hamiltonian 
constraint" 

v 2 

H(x,p) = ^- + V(x)-E^0 . (4) 

Because the canonical Hamiltonian is zero, there are no secondary constraints and 
H is then trivially first class. Jacobi's action in canonical form is thus 

S E [x,p,N}= [ da[px-NH{x,p)] , (5) 

Jo' 

where iV is a Lagrange multiplier. The equations of motion following from variation 
of (5) are 

x = N[x,H} = ^ (6) 
m 

dV 

p = N[p,H] = -N— (7) 

2 

H = — + V-E = . (8) 
2m 
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Combining (6) and (8) determines the multiplier as 

N = x[2(E-V)/m]~ 1/2 . (9) 
The interpretation of (9) is that 

dt = Nda (10) 
is the lapse of physical time, in accordance with the definition of energy. 

The canonical statement of reparameterization invariance is that the action (5) is 
invariant under the gauge transformation given by 

Sx = e[x,H], (11) 
8p = e\p,H], (12) 
8N = e , (13) 

where e(cr') = e(a") = 0. With the choice 

e{a) = (^r^j) T ~ £ daN{a) ' (14) 

where T is the total time 

// 

T= I daN(a) , (15) 

the lapse function is transformed to a constant, namely, N = T/ {a" — a'). This shows 
that every history is gauge related to a history with a constant lapse, and the time T 
is the gauge invariant part of the lapse function. If the histories under consideration 
are restricted to those with constant lapse, the gauge freedom of Jacobi's action is 
removed and (5) becomes 



S E [x,p;T) = I da[px-TH/(a" -a')] . 



(16) 



This form of Jacobi's action is a functional of x(a) and p(cr) and an ordinary func- 
tion of the time interval T. The classical equations of motion for (16), that is, the 
conditions for the extrema of (16), are given by (6) and (7) with iV = T/{a" — a') 
along with 



dS K 1 



°= w = -w^)L daH - (17) 

Since (6) and (7) imply that H is constant, equations (6), (7), and (17) together 
imply Ti. = 0. It follows that the form (16) for Jacobi's action is classically equivalent 
to (5), but has no gauge freedom. 
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4 FUNCTIONAL INTEGRAL FOR JACOBI'S ACTION 

The functional integral associated with Jacobi's action can be constructed by in- 
tegrating over all histories x(a), p(cr), T, with fixed endpoints x(a') = x' and 
x(a") = x" , where the phase for each history is given by the action (16). Thus, 
the functional integral is 



Z E (x", x') = ^— [ dT [ VxVpexpl^ [ da[px - TH/(a" - a')] \ , 

27m J J x (a')=x' I n> J a' J 

(18) 

where VxVp is (formally) the product over a of the Liouville phase space measure 
dx(a)dp(a) / (2ttJi) . The integration measure in (18) can be justified by appealing to 
a BRST analysis based on the canonical action (5), as is done in the Appendix of 
Brown and York (1993b). 

The functional integral over x(a) and p{u) in (18) has the familiar form of the path 
integral associated with Hamilton's action, where a plays the role of time and the 
Hamiltonian is T7i/(a" — a'). This path integral can be written as the matrix 
elements of the evolution operator exp(— iTH/h), so the path integral for Jacobi's 
action becomes 

Z E ( X ", x ') = ^hj dT< x"\e~ lTn/h W > ■ (19) 

Hence, taking the integration of T over all real values, we have 

Z E (x",x') =< x"\6(H)\x' > . (20) 

Note that Z E (x",x') satisfies the time independent Schrodinger equation, namely 
HZ E (x" ', x') = (where H acts on the argument x"), since formally H8(H) = 0. 

From (20) it follows that the trace of Z E (x" , x') yields the density of states 

v(E) = J dx Z E (x, x) = Tr8(H) = Tr5(E - H) , (21) 

where H is the usual Hamiltonian operator. By combining this result with (18) we 
find that v(E) can be written directly as a functional integral, the MCFI: 

v(E) = ^JdTj VxVpexp^ J da [px - TH/2tt] | . (22) 

For later convenience, the endpoint parameter values have been chosen to be a' = —n 
and a" = it. The derivation of the path integral (22) for the density of states shows 
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that the integration can be described as a sum over all phase space curves that begin 
and end at some "base point" x(ir) = x(—tv) = x, plus an integral over the base point 
x. Then roughly speaking, the density of states is given by a sum over all periodic 
histories. However, to be precise, it should be recognized that the sum in (22) counts 
each closed phase space curve a continuous infinity of times because any point on 
the curve can serve as the base point x. Also observe that the integration in (22) 
is over all real values of the time interval T, rather than just positive values. This 
implies that the functional integral for v(E) consists of a sum over pairs of histories 
with members contributing equal and opposite phases (Brown and York 1993b). As 
a consequence, the density of states so constructed is real. 

5 DENSITY OF STATES FOR THE HARMONIC OSCILLATOR 

We now turn to the evaluation of the MCFI (22) for the density of states of a simple 
harmonic oscillator with angular frequency u and Hamiltonian constraint 

9 9 9 

p mu x „ . . 

The periodic nature of the histories suggests the use of Fourier series techniques 
(Feynman and Hibbs 1965) for this calculation. Accordingly, write the phase space 
coordinates as 

oo 

x(a) = ciq + (a,k cos ka + bk sin ka) , (24) 
fc=i 

oo 

p(a) = Co + (c/c cos ka + dk sin ka) . (25) 
k=i 

The functional integral over x(a) and p(a) is replaced by a multiple integral over the 
coefficients in the Fourier series (24) and (25) with measure 

oo 

VxVp = J dao dcQ jfj {dak dbk dck ddk) ■ (26) 
k=i 

Here, J is (formally) the Jacobian of the transformation from x(a), p(a) to ao, cq, 
<2fc 5 bk, Ck, dk- The form (24), (25) of this transformation shows (again, formally) 
that J should be a real constant, and should be independent of T, m, u, and E. (J 
should depend on h, since % appears in the definition of T>xT>p.) One of the goals 
of the present calculation is to determine the real constant J that characterizes the 
change of integration variables specified by (24) and (25). Note that by integrating 
freely over all Fourier coefficients ao, Co, cik, bk, Cf~, dk, we have each closed phase 
space curve correctly counted a continuous infinity of times. This is because the 
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values of the Fourier coefficients depend on the choice of base point that is assigned 
the parameter value a = it (identified with a = —it) on a given closed phase space 
curve. 



With the change of variables (24), (25), the density of states (22) for the harmonic 
oscillator becomes 

v{E) = / dT J da dc (da k db k dc k dd k ) exp{iS E /h} , (27) 

fc=i 

where the phase is obtained by substituting the Fourier series for x(a) and p(a) into 
the action (16): 

<7 — FT m ^ T n 2 T r 2 

- \ Y\^h{a h d k - b k c k ) + £-(4 + d\) + ^-ifll + b 2 k )\. (28) 
fc=i ^ > 

The calculation is simplified by expanding c k and d k about the solutions to their 
"equations of motion". Accordingly, observe that the action (28) is extremized for 
c k and d k that satisfy 

dS T 
= j-^- = ixkb k - —c k , (29) 
oc k 2m 

= || = -nka k - £d k . (30) 
Thus, define new integration variables c k and d k by 

c k = —^-kb k + c k , (31) 
d k = —ka k + d k , (32) 

and the action (28) becomes 

mu 2 T o T 



S„ = ET 



+ dj) + - 4^)(al + bl)\ . (33) 

k=l ^ > 

The integrations over ao, Co, a k , b k , c k , and d k are now straightforward since these 
variables are uncoupled in the action (33). Moreover, for each value of k, the integrals 
over c k are identical to the integrals over d k , and the integrals over a k are identical 
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to the integrals over 6 fe . From these observations it follows that the density of states 
(27) can be written as 



u(E) 



J 



2tvJi 



dT dciodco exp 



mu 2 T 2 T 2 



x 



fc=i 



E[ / da fe dc fe expJ ^ 



m(47r 2 fc 2 - w 2 T 2 ) 2 T _ 2 



4T 



(34) 



Each of these integrals (excluding the integral over T) has the form of a Fresnel 
integral, 

J dxex.p(iAx 2 ) = ^ exp(i7rsigm4/4) , (35) 

where the constant A is real. In evaluating (34) , it is helpful to note that the square 
of the Fresnel integral (35) is iir/A. The result is 



<e) = -u J dT ± n 

Now use the identity 



k=i 



k 



4iv 2 k 2 J 



exp<J -ET 



sin x = x 



fc=i 



n- 



X 



n 



■k 2 



to obtain 



-U 



Next, express the inverse of sin(u>T/2) as 
1 2% 



sin(wT/2) e^ T / 2 - e ~^ T / 2 

oo 

2ie -iuT/2^ e -i«,Tn 



sin(wT/2) 1 [h 



= 2ie~ lujT/2 



exp< —ET 



1 



(36) 



(37) 



(38) 



1 - e-^ T 



(39) 



n=0 

and insert this result into (38). Integrating the series term-by-term, we obtain 

= 2nhJ Y[ ( j ^ 6(E-hu(n+ 1/2)) . 

fc=A " ' n=0 



(40) 



This result shows that the Jacobian J for the change of variables (24), (25) should 
be identified with the (real, infinite) constant 



J 



oo 



2ivh 



k=l 

8 



k_ 

2h 



(41) 



Then the density of states becomes 

oo 

v{E) = ^5{E-hu{n+ 1/2)) , (42) 

n=0 

which is the anticipated result showing that for the harmonic oscillator v(E) is a sum 
of delta functions peaked at half-odd-integer multiples of %uj. 

Finally, we note that the various quantum-statistical and ther mo dynamical proper- 
ties of a system can be obtained from its density of states. In particular, the canonical 
partition function Z{j3) is defined as the Laplace transform of v(E), and from Z{j3) 
the heat capacity, entropy, and other thermodynamical quantities can be found. For 
the harmonic oscillator with density of states (42), the partition function is 

poo 00 

Z{f3) = dE v(E) e~P E = V e -^Hri+i/2) ^ (43) 

^ n=0 

which is the well known result. 
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